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Abstract 

We investigate coherence length effects and hadron attenuation in 
lepton scattering off nuclei in the kinematic regime of the HERMES 
experiment. The elementary electron-nucleon interaction is described 
within the event generator PYTHIA while a full coupled-channel treat- 
ment of the final state interactions is included by means of a BUU 
transport model. The results of our calculations are in good agreement 
with the experimentally measured transparency ratio of incoherent 
electroproduction off ^^N and the multiplicity ratio R\.^ for charged 
hadrons, pions, kaons, protons and anti-protons in deep inelastic scat- 
tering off ^^N and ^^Kr targets. 



1 Introduction 

High energy meson electroproduction off complex nuclei offers a promising 
tool to study the physics of hadron formation. The relatively clean nuclear 
environment of electron induced reactions makes it possible to investigate the 
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timescale of the hadronization process as well as the properties of hadrons 
immediately after their creation. In addition one can vary the energy and 
virtuality of the exchanged photon to examine the phenomenon of color trans- 
parency pn. 

In previous works [21 El El we have developed a method to combine the 
quantum mechanical coherence in the entrance channel of photonuclear re- 
actions with a full coupled channel treatment of the final state interactions 
(FSI) in the framework of a semi-classical transport model. This allows us to 
include a much broader class of final state interactions than usual Glauber 
theory. 

In |21 El we have shown that side feeding effects that are absent in calcula- 
tion based on Glauber theory might play an important role in semi-inclusive 
photoproduction experiments. In ^ we have then investigated the trans- 
parency ratio of exclusive incoherent p° electroproduction off ^^N and ^^Kr 
within our model and found excellent agreement with the Nitrogen data [31 
from the HERMES collaboration. The most important elementary process 
for this reaction was found to be diffractive production on bound nu- 
cleons. Other production mechanisms were supressed by the stringent ex- 
clusivity measure imposed by experiment. We found no signature for color 
transparency in the data for the transparency ratio 

^ " x; 

as a function of the coherence length, meaning that in our model calcu- 
lations we could assume that the diffractively produced p° starts to interact 
with a hadronic cross section right after its production. However, Kopeliovich 
et al. jni have pointed out that the missing signal of color transparency might 
be an accidental consequence of the specific correlation between the of 
the photon and the coherence length of the in the HERMES data. 

A finite formation time of hadrons becomes visible in our investigation 
[T] of the HERMES data on charged particle multiplicities in deep inelastic 
lepton scattering (DIS) off nuclei jHlEl- In the DIS regime we assume that 
the primary production mechanism is governed by the excitation and decay 
of a hadronic string. The fragmentation of the string and the hadronization 
of the fragments takes some time (formation time Tf). After their formation 
time the reaction products undergo hadronic final state interactions with 
nucleons in the target nucleus. Due to time dilatation the formation time tf 
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in the target rest frame can become quite large: 



tj = ^.Tf = Tf. (2) 

Here rrih denotes the hadron mass and is the fraction of the photon energy 
u that is carried away by the hadron. For high energies and for very hght frag- 
ments, e.g. pions, the formation length can easily exceed the size of a nucleus. 
This, however, is not the case for the heavier fragments such as baryons and 
vector mesons for which the formation length at HERMES energies is still of 
the order of the nuclear radius so that "classical" nuclear interactions must 
play a role for these particles. In our model the prehadronic interactions of 
the constituent quarks from the string ends during the formation time are 
effectively accounted for by using the concept of leading hadrons which is 
in accordance with other transport models for high energy reactions such as 
the UrQMD ^] and the HSD model . Alternative explanations for 

the observed modification of the multiplicity spectra in terms of a possible 
rescaling of the fragmentation function in nuclei J3| or a purely partonic 
energy loss ^1 1^ through multiple scattering of the struck quark and in- 
duced gluon bremsstrahlung may model these prehadronic interactions, but 
may neglect the classical FSI after formation. 

The formation time also plays an important role in studies of ultrarel- 
ativistic heavy ion reactions. For example, the observed quenching of high 
transverse momentum hadrons in Au + Au reactions relative to p + p colli- 
sions is often thought to be due to jet quenching in a quark gluon plasma. 
However, the attenuation of high p^ hadrons might also be due to hadronic 
rescattering processes ^Hl if the hadron formation time Tf (in its rest frame) 
is sufficiently short. 

Our paper is structured in the following way: In Sec. 12. II we show how we 
describe the electron-nucleon interaction and how we account for shadowing 
within our model. The transport model is sketched in Sec. 12.21 Our result for 
the transparency ratio of incoherent p° electroproduction is presented in Sec. 
13.11 in comparison with experimental data from the HERMES experiment. 
In Sec. 13.21 we show the calculated multiplicity ratio for charged hadrons, 
pions, kaons, protons and anti-protons. We close with a short summary in 
Sec. H 
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2 Model 



In our approach the lepton-nucleus interaction is sphtted into two parts: 1) In 
the first step the virtual photon is absorbed on a bound nucleon of the target; 
this interaction produces a bunch of particles that in step 2) are propagated 
within the transport model. Coherence length effects in the entrance channel, 
that give rise to nuclear shadowing, are taken into account as described in 
Ref. 0. 

2.1 Elementary photon nucleon interaction and shad- 
owing 

We use the event generator PYTHIA v6.2 ^ITjj to describe the interaction of 
the virtual photon and a nucleon. The basic idea is that in a photon hadron 
collision the physical photon 7 not necessarily interacts as a point particle 70 
(direct interaction) but might fluctuate into a vector meson 
(vector meson dominance, VMD) or perturbatively branch into a qq pair 
before the interaction (generalized vector meson dominance, GVMD). We 
have shown in |4j that the latter is very unlikely in the kinematic regime of the 
HERMES experiment (photon energy v 7 -23 A GeV, pa 0.5-5 GeV^). 
In most cases the photon will therefore interact directly with a parton of 
the target nucleon (DIS) or fluctuate into a vector meson (VMD) before it 
reaches the nucleon. In the latter case the vector meson might either scatter 
diffractively from the nucleon or a hard scattering between the constituents 
of the vector meson and the nucleon might take place. The hard scattering 
like the direct photon nucleon interaction leads to the excitation of one or 
more hadronic strings which finally fragment into hadrons. 

If the struck nucleon is embedded in a nucleus one has to account for its 
Fermi motion and binding energy as well as for Pauli blocking of final state 
nucleons. In addition one has to be aware that the nuclear environment 
influences the VMD part of the photon nucleon interaction, since the vector 
meson components get modified on their way through the nuclear medium to 
the interaction point. The latter has been investigated in detail in Ref. 
There we have shown that at position r inside the nucleus the physical photon 
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state has changed to 
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due to these 'initial state interactions' of the photon. The formfactor Fy 
from Ref. accounts also for contributions of longitudinal photons. If 
one neglects any influence of the FSI for the moment the reaction amplitude 
for the process — > / on a nucleon at position r inside a nucleus changes 
compared to the vacuum in the following way: 



{f\fh) - (/|T|7(r)). 
The nuclear profile function r^\f) can be expressed as 



(4) 



X exp 



-avNi^-iav) / dzkn{b,Zk) 



(5) 



with the nucleon number density n{f), the total vector meson nucleon cross 
section avN and the ratio ay of real to imaginary part of the vector meson 
nucleon forward scattering amplitude. The Bessel function parametrization 
with = 0.78 GeV avoids unphysical contributions from processes where 
Zi ^ Z which would contribute for small values of the coherence length ly- In 
Ref. we have shown that this Bessel function parameterization yields a 
good description of shadowing in photoabsorption down to the onset region 
of shadowing. 

The probability that the photon interacts with a bound nucleon via a 
certain vector meson fluctuation therefore depends on the position inside the 
nucleus. The momentum transfer 



Qv 



— m 



V 
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in the phase factor of Eq. El equals the inverse of the coherence length ly of 
a vector meson component V , i.e. the length that the photon travels as a 
vector meson fluctuation estimated via the uncertainty principle. 
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2.2 Transport model 



The propagation of the produced final state |/) through the nucleus is treated 
within a semi-classical transport model based on the Boltzmann-Uehling- 
Uhlenbeck (BUU) equation. The BUU equation describes the time evolution 
of the phase space density p, t) of particles of type i that can interact 
via binary reactions. Besides the nucleons these particles involve baryonic 
resonances and mesons (vr, 77, p, K, ...) that are produced either in the 
primary reaction or during the FSI. For a particle species i the BUU equation 
can be written as: 

(d dH d dH d \ 
m^Wd-f-Wd-p) mp^^ = icou[fu...f....jM]. (7) 

For baryons the Hamilton function H includes a mean field potential which 
depends on the particle position and momentum. The collision integral on 
the right hand side accounts for the creation and annihilation of particles of 
type z in a collision as well as elastic scattering from one position in phase 
space into another. For fermions Pauli blocking is taken into account in 
Icoii via blocking factors. The BUU equations of each particle species i are 
coupled via the mean field and the collision integral. The resulting system 
of coupled differential-integral equations is solved via a test particle ansatz 
for the phase space density. We also stress that within our model instable 
particles might decay during their propagation through the nucleus. For 
details of the transport model see Ref. (20j . 

The classes of FSI that are included in the transport model go far beyond 
what can be achieved within Glauber theory. As a result the finally observed 
particles do not need to be created in the primary reaction but might be 
produced during the FSI via side feeding. It is therefore clear that a purely 
absorptive treatment of the FSI as in Glauber theory can only be used if one is 
sure that one has eliminated the possibility of side feeding by applying enough 
constraints on the observable (see pj for details). This might be possible for 
exclusive production measured at HERMES but becomes questionable for 
the measurement of inclusive hadron multiplicity spectra. 

As discussed in Sec. I2.1l a hard primary photon nucleon interaction leads 
to the excitation of hadronic strings. The time, that is needed for the frag- 
mentation of the strings and for the hadronization of the fragments, we denote 
as formation time Tf in line with the convention in transport models. For 
simplicity we assume that the formation time is a constant Tf in the rest 
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frame of each hadron and that it does not depend on the particle species. In 
principle, one expects a distribution in the formation times, however, we here 
address only the mean value of such an 'unknown' distribution. We recall, 
that due to time dilatation the formation time t/ in the laboratory frame is 
then proportional to the particle's energy as can be seen from Eq. The 
size of Tf can be estimated by the time that the constituents of the hadrons 
need to travel a distance of a typical hadronic radius (0.5-0.8 fm). We assume 
that hadrons, whose constituent quarks and antiquarks are created from the 
vacuum in the string fragmentation, do not interact with the surrounding 
nuclear medium within their formation time. For the leading hadrons, i.e. 
those involving quarks (antiquarks) from the struck nucleon or the hadronic 
components of the photon, we assume a reduced effective cross section ciead 
during the formation time Tf and the full hadronic cross section ah later 
on. This concept is illustrated in Fig. [TJ due to time dilatation light parti- 
cles emerging from the middle of the string might escape the nucleus without 
further interaction if they carry a large fraction Zh of the photon energy u 
(7 GeV < z/ < 23.4 GeV). However, the hadrons with Zh close to 1 (or energy 
~ u) predominantly stem from the string ends and therefore can interact 
directly after the photon-nucleon interaction. We note, that about 2/3 of 
the intermediate Zh hadrons (mainly pions) are created from the decay of 
vector mesons that have been created in the string fragmentation. Because 
of their higher mass ruh (0.77 - 1.02 GeV) these vector mesons may form (or 
hadronize) inside the nucleus (see Eq. (j2I)) and thus suffer from FSI. The 
effect of the final-state interactions, finally, will depend dominantly on the 
nuclear geometry, i.e. the size of the target nucleus. 



3 Results 

3.1 Incoherent production 

Incoherent p° electroproduction off nuclei has been studied in Ref . |1] . In the 
following we will summarize the results of our investigations. For our model 
calculations we demand that the final state has to consist of two oppositely 
charged pions plus an ensemble of bound nucleons and we use the kinematic 
cuts of the HERMES collaboration This means that we apply the 
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Figure 1: Illustration of an electron nucleus interaction: The virtual photon 
7* excites a hadronic string by hitting a quark q inside a bound nucleon. In 
our example the string between the struck quark q and diquark qq fragments 
due to the creation of two quark-antiquark pairs. One of the antiquarks 
combines with the struck quark to form a 'leading' meson m, one of the 
created quarks combines with the diquark to form a 'leading' baryon h. The 
remaining partons combine to a meson m' that, depending on the mass of 
the meson, might leave the nucleus before it hadronizes (see Q). 

exclusivity measure 

- 2 GeV < = ^1—^ < 0.6 GeV, (8) 
2mN ^ ' 

where rriN denotes the nucleon mass and 

Py = Pn + P-t - Pp (9) 

the 4-momentum of the undetected final state. In Eq. and Pp are the 
4- momenta of the incoming photon and the detected tt+tt^ pair and Pat is 
the 4-momentum of the struck nucleon which, for the calculation of py-, is 
assumed to be at rest. In addition we introduce a lower boundary for the 
four-momentum transfer |t — tmaxl >0.09 GeV^ as imposed by the HERMES 
collaboration to get rid of coherently produced p°. 

In Fig. El we show the transparency ratio Ta for exclusive p° production 
as a function of the coherence length Ip = q~^. The solid line displays the 
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result that one gets if one uses our Glauber expression from Ref. ji] 



a^N^pON d b dzn{b,z) 



X 



1- / dzMb,z,)il-jo{qc\z.-z\))^{l-iap)e"^''^''-'^ 



X exp 



-crpN{l-iap) / dzkn(b,Zk) 



X exp 



-af'' I dz'nib, z') 



(10) 



In Eq. ffTUI) we use for the total p^N cross section ap^ = 25 mb and for 
the elastic part a^' = 3 mb. These two values correspond to the pN cross 
sections used within the transport model for the involved p^ momenta. 

The result of the transport model is represented by the open squares. For 
each data point we have made a separate calculation with the corresponding 
values of u and Q^. In the case of ^^N the Glauber and the transport calcu- 
lation are in perfect agreement with each other and the experimental data. 
This demonstrates that, as we have discussed in Ref. [2], Glauber theory can 
be used for the FSI if the right kinematic constraints are applied. 

After applying all of the above cuts, nearly all of the detected p° stem from 
diffractive p^ production for which the formation time is zero. The ^^N data 
seems to support the assumption that the time needed to put the preformed 
p^ fluctuation on its mass shell and let the wave function evolve to that of a 
physical p° is small for the considered values of Q^. Furthermore, the photon 
energy is too low to yield a large enough 7 factor to make the formation 
length exceed the internucleon distance and make color transparency visible. 
This conclusion is at variance with that reached in Ref. jH]. 

We now turn to ^^Kr where we expect a stronger effect of the FSI. Unfor- 
tunately there is yet no data available to compare with. As can be seen from 
Fig.|21the transport calculation for '^'^Kr gives a slightly smaller transparency 
ratio than the Glauber calculation, especially at low values of the coherence 
length, i.e. small momenta of the produced There are two reasons for 
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Figure 2: Nuclear transparency ratio for electroproduction plotted 
versus the coherence length of the component of the photon. The data 
is taken from [S]. The solid line represents the Glauber result when using 
Eq. (fTn|l . For each transparency ratio calculated within our transport model 
(open squares) we used the average value of and v of the corresponding 
data point. The experimental exclusivity measure (jHI) has been taken into 
account for both nuclei. 

this: About 10% of the difference arises from the fact that within the trans- 
port model the p° is allowed to decay into two pions. The probability that at 
least one of the pions interacts on its way out of the nucleus is about twice 
as large as that of the p*^. The other reason is that in the Glauber calcula- 
tion ()10|1 only the inelastic part of the p^N cross section enters whereas the 
transport calculation contains the elastic part as well. Thus all elastic scat- 
tering events out of the experimentally imposed t-window are neglected in 
the Glauber description. It is because of this t-window that also elastic p^N 
scattering reduces the transport transparency ratio shown in Fig. |21 Both 
effects are more enhanced at lower energies and become negligible for the 
much smaller '^^N nucleus. 
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3.2 Hadron formation and attenuation in DIS 

As mentioned above hadron production in deep inelastic lepton-nucleus scat- 
tering (DIS) offers a way to study the physics of hadronization The 
reaction of the exchanged virtual photon (energy u, virtuality Q"^) with a 
bound nucleon leads to the production of several hadrons. Depending on the 
space-time evolution of the production process the rescattering in the sur- 
rounding nuclear medium will change the energy distribution and the mul- 
tiplicity of the produced particles. Consequently, the particle spectrum of a 
lepton-nucleus interaction will differ from that of a reaction on a free nucleon. 
In order to explore such attenuation effects the HERMES collaboration has 
investigated the energy u and fractional energy = Eh/u dependence of the 
charged hadron multiplicity ratio 



in DIS off N jH] and Kr ^ nuclei. Here Nh{z^v) represents the number of 
semi-inclusive hadrons in a given (z, z/)-bin and Ne{v) the number of inclusive 
DIS leptons in the same z/-bin. 

It was suggested in Ref. jH], that a phenomenological description of the 
i?5vf data can be achieved if the formation time, i.e. the time that elapses 
from the moment when the photon strikes the nucleon until the reaction 
products have evolved to physical hadrons, is assumed to be proportional to 
(1 — z/i)z/ in the target rest frame. This {1 — Zh)v dependence of the formation 
time Tf is compatible with the gluon-bremsstrahlung model of Ref. |24]. In 
the investigations of Ref. [Hj any interaction of the reaction products with 
the remaining nucleus during this formation time has been neglected. After 
the formation time the hadrons could get absorbed according to their full 
hadronic cross section. Another interpretation of the observed BJIj spectra - 
as being due to a combined effect of a rescaling of the quark fragmentation 
function in nuclei due to partial deconfinement as well as the absorption of 
the produced hadrons - has recently been given by the authors of Ref. [13] . 
Furthermore, calculations based on a pQCD parton model IT^ explain 
the attenuation observed in the multiplicity ratio solely by partonic multiple 
scattering and induced gluon radiation. It has already been pointed out by 
the authors of Ref. JH] that a shortcoming of the existing models is the 
purely absorptive treatment of the final state interactions (FSI). We avoid 
this problem by using the coupled-channel transport model. 




(11) 
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In our calculation we employ the kinematic cuts of the HERMES exper- 
iment as well as the geometrical cuts of the detector. In actual numbers: we 
require for the Bjorken scaling variable x = 2mi^u ^ ^-^^ (with niN denoting 
the nucleon mass) , for the photon virtuality > 1 GeV^ and for the energy 
fraction of the virtual photon y = u/Efyeam < 0.85. In addition, the PYTHIA 
model introduces a lower cut in the invariant mass of the photon-nucleon sys- 
tem at W = 4 GeV that is above the experimental constraint W > 2 GeV. 
This limits our calculations to minimal photon energies of Umin = 8.6 GeV 
as compared to Umin = 7 GeV in the HERMES experiment and leads to a 
suppression of high events at energies below u ^ 15 GeV. 

We have shown in Ref. [7] that our model simulation reproduces the 
experimental average values of the kinematic variables and u {zh) on A'^and 
Kr as a function of Zh {i') very well. Since the particles with Zh close to 1 are 
predominantly leading hadrons we could use the high Zh part in the fractional 
energy spectrum to fix the leading hadron cross section. A good agreement 
with the data is achieved for aiead = 0.33ah during the formation time Tf. 
We note, that this value for aiead represents an average value over time from 
the virtual photon-nucleon interaction to the actual hadron formation time. 
Since for light nuclei only a fraction of the leading hadrons is formed inside 
the nucleus, the effective leading hadron cross section aiead is expected to 
be smaller accordingly. For a detailed investigation we refer the reader to a 
forthcoming study j2Sl- 

In Fig. El we show the influence of different formation times Tf on i?^ 
using the effective cross section aiead = 0.330"^. We find, that formation 
times Tf > 0.3 fm/c are needed to describe the experimental data with 
little sensitivity to higher values. This is compatible to the range of values 
extracted from the antiproton attenuation studies in Ref. |2tij . 

Some of our model assumptions, e.g. the local density approximation, 
become questionable for very light nuclei. We have, therefore, used the Kr 
data to fix the value of the formation time. However, we also get a satisfying 
agreement with the Nitrogen data. 

The pQCD model of Ref. ^3] predicts a hadron attenuation ~ A^^^ since 
the parton energy loss is proportional to the propagation length squared. It 
is thus important to get the scaling with target mass A from our present 
approach in order to allow experimental studies to distinguish between the 
different concepts. To this aim Fig |21 also shows predictions for a Xe target. 
In accordance with the authors of Ref. we predict only a small change 
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Figure 3: Calculated multiplicity ratios of charged hadrons for A^, Kr and Xe 
targets for a fixed leading hadron cross section aiead = O-SSa^ and different 
values of the formation time Tf. The data for the Nitrogen target have been 
taken from Ref. jH] while the preliminary Kr data stem from Ref. [H]. 

in the multiplicity spectra compared to the Kr target such that the scaling 
exponent is lower than 2/3. 

In Fig. E]we show the results for the calculated multiplicity ratio of tt~ , 
7r+, K~ , K~^, p and p for Kr in comparison with the experimental data. In 
our calculations we use again the kinematic cuts of the HERMES experiment 
with the additional restrictions x > 0.02, E.,^j^=2.5-15 GeV and £"^,^=4-15 
GeV that are experimentally necessary for particle identification j2Z|. We 
use a constant formation time of 0.5 fm/c and aiead = 0.33(j/j for all hadrons. 
Without further fine tuning we get a satisfying description of all the data 
meaning that the formation times of mesons, baryons and antibaryons are 
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Figure 4: Calculated multiplicity ratios of vr"*", vr", , K~ , p and p for 
Krypton using a fixed leading liadron cross section aiead = 0.33ah and for- 
mation time Tf = 0.5 fm/c. The experimental data has been taken from Ref. 

m- 

about equal. 

4 Summary 

We have developed a method to account for coherence length effects within 
a semi-classical BUU transport model. This allows us to describe incoherent 
meson photo- and electroproduction off nuclei at GeV energies using a full 
coupled channel treatment of the FSI. We have calculated the transparency 
ratio for exclusive incoherent p° photoproduction off ^^N and ^^Kr. The re- 
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suit for ^^N is in agreement with experimental data and with the Glauber 
prediction. The latter shows that in the case of ^^N Glauber theory is appli- 
cable after the kinematic cuts of the HERMES experiment are applied. Since 
we do not use a formation time for diffractively produced vector mesons we 
deduce that one cannot see an onset of color transparency in the Nitrogen 
data. For the '^^Kr target no experimental data is available to compare with. 
However, we find deviations from the simple Glauber model because of the 
finite life time of the p° and elastic scattering out of the kinematically allowed 
|t|-region. As discussed by Kopeliovich et al. jH] one might see an onset of 
color transparency when investigating the transparency ratio as a function 
of for fixed coherence length. 

In addition we have shown that one can describe the experimental data 
of the HERMES collaboration for hadron attenuation on nuclei without in- 
voking any changes in the fragmentation function due to gluon radiation. 
In our dynamical studies, that include the most relevant FSI, we employ 
only the 'free' fragmentation function on a nucleon and attribute the hadron 
attenuation to the deceleration of the produced hadrons due to FSI in the 
surrounding medium. We find that in particular the ^-dependence of R\j is 
very sensitive to the interaction cross section of leading hadrons and can be 
used to determine aiead- The interaction of the leading hadrons during the 
formation time could be interpreted as an in medium change of the fragmen- 
tation function. The extracted average hadron formation times of Tf > 0.3 
fm/c are compatible with the analysis of antiproton attenuation m p + A 
reactions at AGS energies [21]. 

In this work we have also for the first time compared our results for the 
attenuation of vr^, 7r+, , K~^, p and p with the experimental Kr data. We 
find a satisfying description of all data using the same constant formation 
time Tf=0.5 fm/c and the same reduction of the effective leading hadron 
cross section cT/eacZ=0.33 fm/c for all hadron species. 
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